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• Abstract 

£h ' We observe chaotic dynamics in a damped linear oscillator, which is 

driven only at certain regions of phase space. Both deterministic and 
random drives are studied. The dynamics is characterized using standard 
techniques of nonlinear dynamics. Interchanging roles of determinism and 
^v^j , stochasticity is also considered. 
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CM ! 1 Introduction 

00 

C*") ' It is generally known that continuous linear systems cannot exhibit chaotic 

. time evolution. However recently it has been shown that chaotic behavior can 

occur in systems without explicit non linearities. Hirata et.al. [I] have shown 
that continuous chaotic wave forms can be constructed by the superposition of 
certain pulse basis functions. In a linear second order filter driven by randomly 
polarized pulses, chaotic behavior is observed when the time series is viewed 
backwards in time [2]. Further they have shown that by reversing the time 
series, folded band chaos similar to Rossler attractor can also be synthesized 
[3] . Even if the time reversal is not physical, this reveals the importance of pure 
randomness which is associated with chaotic behavior. Studies in this direction 
are important in understanding the relation between randomness, determinism 
and complexity in dynamical systems. 

We study a weakly damped driven linear oscillator where the oscillator and 
the drive are coupled only when the trajectory is in a thin strip in the phase 
space. This modification introduces strong nonlinearities required for chaotic 
behavior. First we consider the oscillator coupled to a deterministic drive. Then 
we consider a random forcing. We show that the attractors formed by both the 
deterministic and random forcing are similar in many respects. 
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2 Basic model 



The driven damped linear oscillator is well known and is important in almost all 
branches of physical science. The dynamics of this system satisfies the equation, 

^-+"/^ + fl 2 x = asm(u>t). (1) 
dt z at 

The parameters in the equation are 7, the damping coefficient, f2, the angular 
frequency of the oscillator and to the frequency of the drive. Such a system 
shows periodic behavior after the initial transients have vanished. 

Our basic model is the one in which the system is driven only in a thin strip 
of the phase space. We choose the origin as the location of the strip because 
orbits with lower energy may not encounter the strip, if it is located far from 
the origin. Such orbits will spiral towards zero as in the case of an ordinary 
damped oscillator. When not driven, the oscillator performs a damped motion 
and the drive runs freely. First we consider a deterministic drive. The dynamics 
is characterized using standard techniques. Also the behavior under parameter 
variations is studied. Later we discuss the effect of replacing the deterministic 
driving term with a stochastic one. 



3 Deterministic driving 

3.1 Dynamical equations 

The equation [T] modified accordingly by the addition of a new term n(s,x) as, 

+ 7""77 + ^ 2 ' x = fJ-(s,x)asm(wt) (2) 
dt z at 

where, 

I \ - / 1 if I x I - s d\ 

- I o otherw ise. [6 > 

Thus with a deterministic drive, what the dynamical equations represent is 
legitimate deterministic system in every sense. 



3.2 Numerical results 

Numerical simulations were done using the Runge-Kutta algorithm with a step 
size of 10~ 3 . We choose the following values for the parameters: f2 = 1.0, 
7 = 0.1, to = 0.3, s = 0.1 and a = 2.5. The parameters were chosen 
by numerical trials to express the concepts clearly, and are fixed during the 
evolution of the system in time. Fig[T]shows the phase space plot of the system. 
Note that there is a discontinuity in a strip of width 2s near in the the phase 
space, where the oscillator is coupled to the drive. A closer view of the strip is 
given in Fig. [4] It can be seen that the trajectories that approach the origin 
may either cross the strip with or without considerable modifications, or it may 
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get reversed. Fig. [5] shows the time series x(t) of the system. The amplitude of 
oscillations, and the number of 

oscillations between two consecutive crossings of the strip are irregular. This 
corresponds to the phase space dynamics which consists of reflections in the strip 
and transitions which modify the trajectories. The power spectrum as shown 
in Fig. [3] is broad which confirms the aperiodic behavior in time. Varying the 
parameters 7 and s, the system exhibits both chaotic and periodic oscillations. 
From the bifurcation diagrams shown in Fig. [5] and Fig. it can be seen that 
the chaotic behavior disappears for higher values of 7 and s. Also note that for 
s ~ 1 there is a transition to the case of an ordinary driven damped linear 
oscillator, where no chaotic behavior is expected. 

We used DATAPLORE [4] for the analysis of the time series. The time de- 
lay required for reconstructing a timeseries is a widely discussed topic [5] . With 
a delay of 1.56, it is found that the autocorrelation function falls to 0.5. With 
this delay, embedding dimension was calculated using False nearest neighbors 
method 6J. It is found that the percentage of nearest neighbors approaches 
zero with embedding dimension 3. Lyapunov Exponents were calculated with 
embedding dimension 3, number of nearest neighbors 37 and degree of the ex- 
trapolating polynomial 3. The obtained result is multiplied by the sampling 
frequency of the input signal to calculate the Lyapunov exponents. The Largest 
Lyapunov exponent is found to be 3.19. We obtained 2.148 as the Kaplan- Yorke 
dimension. The reconstructed phase space of the attractor is shown in fig. [7] 
and is similar in appearance to other chaotic systems. 

4 Random driving 

4.1 Dynamical equations 

Random forcing is achieved by the drive assuming a random value £ each time 
the phase space trajectory enters the strip. The amplitude of forcing £ is Gaus- 
sian random variable with variance a 2 = A- and zero mean. We define x(s, x) 
as follows, 

X[ ' 1 \ otherwise. [ ' 

Retaining fj,(s,x) defined as in eq. [5] for comparison, the dynamical equation is, 

+ 1 ~<E + ^ 2x = a f i ( s ^ x )^(s,x) (5) 

4.2 Numerical results 

When the oscillator is driven randomly, it is found that the attractor is still 
chaotic with the largest Lyapunov exponent equal to 2.92 and a Kaplan- Yorke 
dimension of 2.31 when embedded in three dimensions. The dynamics in phase 
space is similar to the one obtained with a deterministic drive as shown in Fig. 
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[5J The reconstructed attractor (with time delay At = 1.1) given in Fig. [5] is 
similar to the one obtained with a deterministic evolution. 

5 Discussion 

The effect of a forcing on the dynamics can depend on many factors. The most 
important are the amplitude and the time for which such a forcing persists. 
With Eq. [2] the dynamics of the system in terms of phase space variables is 
well defined. The only degree of freedom available for an irregular behavior 
is through the time interval in which the system is driven and is not driven. 
At the onset of chaos, it can be seen from Fig. [TO] that the durations for 
which the driving persists are highly irregular. Also irregular are the intervals 
in which the system is not driven. This intermittent nature of the drive and 
irregular durations of the driving events are the only source of irregularities 
for a dynamics which is linear in every region of phase space. From Fig. 1111 
it can be seen that the random forcing applied to the oscillator also persists 
for irregular durations. Stochasticity as such is a source for irregularities in a 
system. But it may not be the same as the topological structure of an attractor 
demands. There, the intermittent driving and irregular driving durations play a 
significant role. It modifies the irregularities associated with the driving suitable 
for chaotic dynamics in the system. The durations for which the system is not 
driven are also irregular in a similar manner. Thus intermittency is significant to 
the chaotic behavior due to its role in synthesizing or seasoning the complexities 
that occur in the dynamics. 

The damped oscillator, together with its drive, have got another very funda- 
mental similarity with deterministic chaos. The phase space of a chaotic system 
is dense in periodic orbits, or, it contains infinitely many periodic orbits. None 
of these periodic orbits are stable, but they significantly influence the evolution 
of a system in the phase space. This aspect is important in considering the 
anatomy of the chaotic behavior exhibited by the system under consideration. 
Different periodicities arise due to different numbers of reflections and trans- 
missions in the strip. The drive acts as the triggering mechanism for reflections 
or transmissions in the strip which facilitates smooth transitions between or- 
bits of various periodicities. This is also the reason why the attractor structure 
remains similar under both deterministic and periodic drivings. Though phys- 
ically similar, two driving schemes are completely different in a mathematical 
sense. Thus intermittently driven oscillator is an example where determinism 
can mimic stochasticity. 

6 Conclusion 

A Damped linear driven oscillator can exhibit chaotic behavior if the coupling 
between the drive and the oscillator are applied only at a narrow region of the 
phase space. Similar chaotic behavior can be observed with a purely random 
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forcing. Though physically similar mechanisms exist for such a behavior, math- 
ematically, they are different. We hope that this work has revealed yet another 
relation between randomness and determinism in nature. 
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Figure 3: Power spectrum of the time series x(t) (f2 = 1.0,7 = 0.1, lo = 0.3, 
,s = 0.1 and a = 2.5). 
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Figure 4: Close view of the strip. Note the trajectories: a) reverse direction b) 
cross the strip but suffers deformation, c) cross the strip without considerable 
modification. 
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Figure 5: Bifurcation diagram obtained by varying the strip width s. Chaotic 
behavior disappears as the width of the strip is (Q = 1.0, 7 = 0.1, u = 0.3 
and a = 2.5). 




Figure 6: Bifurcation diagram obtained by varying the damping coefficient 7. 
Chaotic behavior is suppressed as the damping is increased (ft = 1.0, u = 0.3, 
s = 0.1 and a = 2.5). 
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Figure 7: The delay reconstructed attractor (ft = 1.0, 7 = 0.1, w = 0.3, 
s = 0.1, a = 2.5 and At = 1.56 ). 




Figure 8: Close view of the phase space strip of a randomly driven oscillator. 
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Figure 9: The reconstructed attractor resulting from random driving (ft = 1.0, 
7 = 0.1, s = 0.1, a = 2.5 and At = 1.1). 
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Figure 10: Deterministic drive: the duration of the driving episodes. N repre- 
sents N th driving event. 
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Figure 11: Random forcing: the duration of the driving episodes. N represents 
N th driving event (At = 0.25s, Q = 1.0,7 = 0.1, s = 0.1 and a = 2.5). 
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